Kondo Effect in Fermi Systems with a Gap: A Renormalization Group Study 
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We present the results of a Wilson Renormalization Group study of the single-impurity Kondo 
and Anderson models in a system with a gap in the conduction electron spectrum. The behavior of 
the impurity susceptibility and the zero-frequency response function, T « S z ; S z » axe discussed 
in the cases with and without particle-hole symmetry. In addition, for the asymmetric Anderson 
model the correlation functions, < S ■ 5(0) >, < >, and < na(2 — n^) > are computed. 
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I. INTRODUCTION 



The properties of a magnetic impurity in a semiconductor or an insulator are of interest for a variety of reasons. In 
a normal fermi system a spin-i impurity yields logarithmic temperature dependences in the impurity susceptibility 
and the resistivity at high temperatures; at low temperatures the magnetic moment is quenchedElu. The existence 
of a sharp fermi surface and the concomitant occurrence of (low-energy) particle-hole pairs play an important role 
in understanding the behavior of the model. Thus it is interesting from a theoretical point of view to understand 
whether the Kondo effect persists and under what conditions quenching occurs in a system with a gap and determine 
the behavior quantitatively. We also note that the Anderson impurity has been studied in the context of the loga- 
rithmic temperature dependence of the conductivity of trans-polyacetyleneH; the system was modeled by a continuum 
Hamiltonian that exhibits a gap due to Peierls distortion. The impurity model was investigated using a Hartree-Fock 
closure of the equation of motion. In addition, a variety of Kondo and valence fluctuating insulators (modeled theo- 
retically by a Kondo or Anderson lattice) such as SitiBq and Ce^Bi^Pt^, among others provide another motivation 
I ' for studying the single impurity problem in a system with a gap. 
, In this paper we present the results of our study of the Kondo and Anderson impurities in a system with a gap. We 
■ apply Wilson's (numerical) Renormalization Group (RG) technique using a variant of a numerical tridiagonalization 
t**"* ' method devised by us earlienj and provide results for both the susceptibility and zero-frequency response functions. 
We also discuss a simple effective Hamiltonian that allows us to understand the physics underlying our results. 

We begin with a summary of previous work on the problem; in the next section we provide a brief sketch of the 
technique that is described in detail in the literatures. We then present the results of our numerical simulations and 



' in the final section we discuss the effective Hamiltonian description of our results. UTl 

We begin with, a . b rief overview of previous work. The first calculations were done by Ogura and Sasouilj and 
Takegahara et alBM. O gura and Saso used a 1/JV expansion of the degenerate Anderson model and found to leading 
order a transition between the triplet and singlet ground states when the gap E g equals twice the Kondo temperature 
I ■ Tjf. In their Quantum Monte Carlo (QMC) simulations they found indications of a similar transition even for the 
symmetric Anderson modela; for the asymmetric model they obtained a transition between the different ground states 
at approximately E g s» 3Tk- We note that their QMC computations were limited to temperatures above Tk/10. 

Takegahara, Shimizu, and SakaiErEj used both Quantum Monte Carlo simulations and the numerical Renormaliza- 
tion group method of Wilson. They considered the symmetric Anderson model and found that at low temperatures 
the susceptibility follows a Curie law resulting from an unquenched magnetic moment. They observed the crucial 
difference between symmetric and asymmetric Anderson models cases; when particle-hole symmetry is obeyed the 
moment remains unquenched for all non-zero values of the gap while there is a transition in the asymmetric case. 
They also used the Wilson numerical Renormalization group to follow the spectrum of the low-energy states but not 
the susceptibility. It is difficult to use their version of the numerical RG formulation to calculate low temperature 
(much less than thp-band gap) properties of the model. 

Yu and Guerrerdia studied an Anderson impurity in a semiconducting host using the density matrix renormalization 
technique. Their calculation which is restricted to T — considered electron spin-impurity spin correlation functions 
and found no qualitative difference between the symmetric and asymmetric cases. We will .comment on this point 
later. The importance of the particle-hole symmetry breaking has been emphasized recentlyt£l in the context of the 
Kondo problem with a pseudogap. In this work the impurity susceptibility for the case of a Kondo system with a gap 
was also calculated using the Wilson renormalization group method: in the particle-hole symmetric case the impurity 
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retains its moment in the ground state for all </; in the presence of potential scattering the moment is completely 
quenched provided that A << Tk- These results are in agreement with those of Takegahara et a£3 ! E3. 

We have made a comprehensive study of the Kondo and Anderson models in gapped systems with and without 
particle-hole symmetry breaking using the .Humerical RG method. Our RG formulation is based on the numerical 
tridiagonalization technique developed by usta, which allows us to calculate various quantities in the entire temperature 
range. We report results for a zero-frequency response function, correlation functions and the susceptibility: we 
emphasize the differences in their behaviors in the various regimes and clarify which of these are good probes of the 
nature of the low-temperature fixed point behavior. 



II. WILSON'S RG FORMULATION 



We consider the Kondo and Anderson models with a conduction electron Hamiltonian with the density of states 
p(e); as a function of the energy e, p is a constant for Do > \e\ > Aq, where the band edges lie at ±Do from the Fermi 
level which is chosen to be in the middle of the gap. The width of the gap is thus 2Ao. The impurity part of the 
Hamiltonian for the spin— i Kondo problem with the impurity spin denoted by S is given in standard notation by 
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where cj~ creates an electron with spin p at the impurity placed at the origin. The choice K = in the Kondo 
problem and + U- = in the Anderson model correspond to particle-hole symmetry. 



Following Wilson we perform a logarithmic discretization of the energy variable; we rescale the energy by Dq so 
that e € [—1, 1], introduce a scale factor A(> 1), and define the nth interval for positive e to lie between A - " -1 and 
A~™. The band gap is chosen to be 

A = A~ Mo . 

Next we replace the continuous set of energy levels in the nth interval [A - ™ -1 , A~™] and [— A~™, — A - ™ -1 ] by 
single levels at (A _n_1 +A~™)/2 and — (A~™+A~™ -1 )/2 respectively, and introduce a+ M , 6+ , the conduction electron 
creation operators for the states with the corresponding energies (A _,l_1 +A _n \ 7 /2 and — (A - ™ -1 + hr n )/2. After this 
discretization, the Anderson Hamiltonian can be rewritten in the following formd (Here we only consider the Anderson 
Hamiltonian; similar RG formulation can be written down for the Kondo Hamiltonian) 
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The initial values of the couplings T (= 0.57rl^ 2 ), e^, and U, are now in units of Dq (taken to be one); the gap in the 
density of states is between A and — A. p, 

We use the (numerical) tridiagonalization scheme devised by us earlien3 to transform the Hamiltonian to the 
following tridiagonalized form: 



1 + A 



-1 No-l 



n=0 



+ 



2T 

7T 



(4) 



where iV = 2M - 1. 

In order to carry out the RG calculation, we need to rescale the Hamiltonian at each iteration step. The rescaling 
is done by defining Hjy as follows: 
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where e 
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The recursion relation can be written in the following compact form: 
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This recursion relation enables one to set up an iterative diagonalization scheme to calculate the energy levels of Hn 
and thus to determine thermodynamic properties; the recursion is implemented numerically and is stopped at N = Nq 
corresponding to the edge of the gap below which there are no conduction electron states. Recall that as we increase 
N, the system effectively evolves from high temperatures to low temperatures. At a given N, the thermodynamic 
quantities are calculated for Tn = 1/(/3Sn) for a selected values of /?. By studying the evolution of the many-body 
energy level structures we also obtain information near the fixed points of the Hamiltonian. 

For N < N , the thermodynamic quantities are calculated for = l/(/3SW) for a selected value of 0; the 
accuracy of the numerical evaluations is enhanced by performing a second-order perturbation calculation by writing 
the Hamiltonian as 



where 



and 



Hp 



H A = (H N + H I + H B )/S N , 



Hi = -^—{ftrJN+iti + h.c), 



1 f N °~ 1 



For N = Nq, the thermodynamic quantities are calculated for a sequence of temperatures {T;}. Since H^ is the 
full Hamiltonian (hence, no second order perturbation is needed), we can calculate the quantities at temperatures 
much lower than typical energy scale at N = Nq, which is the bandgap A. We choose 7] to be a sequence of values 
from 0.175 of the maximum energy kept in the many body states of Hnq to 0.000175 of the maximum energy. Thus 
the thermodynamic quantities at low temperatures are calculated with the "effective Hamiltonian" Hn ■ 



III. RESULTS 



We present the results obtained from our numerical calculations for the two models. 
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A. Kondo Model 



Our calculations were performed using a scale factor of A = 2 and a band gap energy A = 1.22 x 10~ 4 corresponding 
to Mo = 13 . The first and obvious quantity to consider is the impurity susceptibility, x, which we emphasize is defined 
as the total susceptibility of the system minus the susceptibility of the pure system. In the Kondo problem in the 
absence of potential scattering T\ approaches the value 1/4 as T — > for any finite bandgap. The ground state 
is a magnetic doublet, its quantum numbers are (Q — 0, S — 1/2). This is in agreement with the results of 
Takegahara et all3 for the symmetric Anderson model. The susceptibility curves are displayed in Fig. 0(a). Note 
that some data obtained at intermediate points have been suppressed for clarity in this figure as well as in other 
figures we are going to present in this paper. The calculation is done for initial values of the coupling given by 
Jo = —0.1, —0.2, —0.3, —0.4, —1.0. Note that for large values of \ J\ the universal shape of T\ of the ordinary Kondo 
problem is evident at high temperatures but below temperatures of the order of the gap T\ increases sharply. 

The, effect of particle- hole symmetry breaking introduced by potential scattering is very important as has been noted 
beforefcJi 2 ]. The results for T X are displayed in Fig. |(b). For K = 0.1 and J = -0.2 (T K « 7.4 x 10~ 6 << A), T X 
again goes to 1/4 as T goes to zero. For stronger Kondo coupling, J = —0.4 (Tk ~ 2.1 x 10~ 3 >> A), the impurity 
spin is quenched and T\ — > 0. There is a discontinuous ("first-order") transition due to a crossing of energy levels. 
Crudely speaking, in the generic case without particle-hole symmetry, the transition occurs when the energy gained 
by forming the singlet which is of the order of Tk is larger than the energy required to create a particle-hole excitation 
across the gap. 

We have also calculated the zero- frequency response function T « S Z ;S Z >>. The techniques for performing 
such calculations have been explained in an earlier paper for the ordinary Kondo problemll We used A = 3.0 in the 
calculation of the response function. Note that for the ordinary Kondo problem << S z ; S z >> is essentially the same 
as the impurity susceptibility x for small values of the initial coupling!. However, for the density of states with a 
gap, << S z ; S z » and x behave quite differently at low temperatures (when T < A). In the absence of potential 
scattering, in contrast to T\ which approaches a fixed value of 1/4 as T — > 0, T « S z ; S z » approaches a value 
Co which depends on the bandgap; this persists also when the moment is not quenched in the presence of potential 
scattering. The results for Co are listed in Table I. For K = our results are consistent with the value Co being 
proportional to A 2 for A << Tk- We will derive this result from our effective Hamiltopian description in the next 
section. This result agrees with the claim made by Takegahara et al. for the susceptibility^; we note that they appear 
to have identified << S Z ;S Z >> with the impurity susceptibility. In the presence of potential scattering when A is 
increased for fixed Tk the ground state changes abruptly from a singlet [(Q = —1, S — 0)] to a doublet. The value of 
T\ jumps from to 1/4 and correspondingly the value of Co also jumps discontinuously. 

B. Anderson Model 

The calculations for the Anderson model were performed with the parameter A = 3. A range of values was used 
for the band gap A = A _M °: the value of Mq was varied between 3 and 19. 

For the symmetric Anderson Model, with U = 0.1, e<j = — ^ = —0.05, and T — 0.006, T\ reaches the value of 1/4 
as zero temperature is approached irrespective of the value of the bandgap A, signaling a doublet ground state and 
an unqucnchcd impurity moment. If A << Tk ~ 5.12 x 10 -6 , T\ first decreases toward zero along the universal 
Kondo curve; however, when T < A, it rises to 1/4 as T goes to zero. If A is comparable or larger than Tk, on the 
other hand, T\ gradually increases to 1/4. Our results for T\ are displayed in Fig. |](a). 

The case of the the asymmetric Anderson Model was studied using the parameter values U = 0.1, = —0.0001, 
and T = 0.00015, and the results are displayed in Fig. ^(b). When A = 0, the system goes successively-through 
the free-orbital regime, the mixed valence regime, the local moment regime, and the frozen moment regimea. When 
A 7^ 0, for T > A, the T\ curve initially follows the curve for A = as the temperature is lowered. When T drops 
below A, T\ curves starts to rise. For A > Tk, the curve continues to rise to 1/4 as T goes to zero. On the other 
hand, when A < Tk, the curve stops rising, turns over and tends to zero as T — ► 0. This behavior is clearly similar 
to that of the K ^ case of the Kondo model. 

Fig. U shows the general temperature dependences of the zero-frequency response function T « S z ; S z » for the 
asymmetric case. There is no qualitative difference in the behavior of the response function T << S z ; S z » between 
the Kondo and Anderson models. In the symmetric case where the ground state is characterized by (Q = 0, S — 1/2), 
we again found that the zero-temperature value Co is proportional to A 2 when A decreases. Also for the asymmetric 
case Co jumps discontinuously as the ground state changes from a singlet to a doublet as A is increased. 

In addition, we have also computed the following correlation functions: < S ■ <?(0) >, < >, and < n^(2 — n^) >■ 
Representative figures are shown in Fig. ||(a), and Fig. ||(b). Here the main point to be emphasized is that once there 
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are no charge fluctuations (for example, when the system approaches the local moment regime, or when T < A) the 
correlation functions do not change and approach constant values. In particular, when the local moment regime is 
reached (A is less than the temperature for the local moment formation), the correlation functions tend to the same 
constants as T — > independent of the band gap. While the mixed valence regime is still reflected in the temperature 
dependence of the correlation functions, the Kondo effect does not show up in the correlation functions. This point 
is not very well appreciated. One simply cannot investigate the Kondo effect using local correlation functions, such 
as the impurity spin-conduction electron spin-density at the origin, since they do not contain information about the 
system on the energy scale of Tk or equivalently the length scale of hv e /TK where v e is.-the characteristic velocity 
of the electrons. We believe that this is the reason why in the work of Yu and Guerrerdiil no difference was found 
between the symmetric and asymmetric Anderson models in the correlation functions at short length scales. 

Finally, we present our results for the mixed-valent regime. We considered the asymmetric Anderson Model, with 
U = 0.1, €d = —0.025, and T = 0.01. When A = 0, the system goes from the free-orbital regime through the mixed 
valence regime directly to the frozen moment regime, without going through the local moment regime. The results 
for the susceptibility, < >, and T « S z ; S z » are shown in Fig. ||(a)-(c). Again, depending on the value of the 
band gap, T\ can go to zero or 1/4 (there is a sharp transition). For the cases that T\ goes to zero, T « S z \ S z >> 
also goes to zero, and all correlation functions approach constants, which are independent of the band gap. But for 
the cases that T\ goes to 1/4, both T « S Z ;S Z >> and the correlation functions approach the values which are 
band-gap dependent. 

IV. EFFECTIVE HAMILTONIAN DESCRIPTION 

In this section we provide a simple interpretation of the low-temperature behavior of the models in the various 
regimes on the basis of a simple effective Hamiltonian. Let us consider first the Kondo Model with the gap in the 
conduction electron density of states between — A to A. The initial couplings are Jo and Kq in units of the bandwidth 
Do, which is taken to be unity. In our RG calculation, the band gap is taken to be A = A~ M ° , where Mq is an integer; 
this corresponds to the maximum N being iVo = 2Mo — 1 — there are even number of conduction electron levels in 
the discretized system. 

Imagine that we have successively integrated out the high energy degrees of freedom and arrived at the effective 
Hamiltonian at the energy scale A = A~( No+1 ^ 2 ; as we pointed out earlier the iterative RG procedure cannot be 
carried beyond this energy scale corresponding to the maximum iteration number No since there are no conduction 
electron states left. The low temperature properties (i.e, for T << A) can be calculated with this effective Hamiltonian. 

Let us consider the case when A >> Tk- The effective Hamiltonian is close to that of the J — fixed point and 
can be written, keeping the leading order terms, as 

H eff = -JS ■ a(0) + Kf+f^ + A(a+a^ - b+b^) . 

Here <?(0) = ^f^^nufv and = -^( a ti + fyt)- hi the above effective Hamiltonian we have only kept the lowest single 
electron/hole levels of the conduction electron Hamiltonian; these are represented by the creation operators a+ and 
b+ and we have neglected the irrelevant operators. Since fo is proportional to A~ N °/ 4 , we have§ 

fa =a Q A- No/4 (a^ + b^ + .... 

Thus the first two terms of H e ff are marginal, and J and K must scale as J = JqA and K — KoA. When | Jq\ << 1 
and |-Ko| << 1, the last term dominates, and the ground state is a doublet. 

Next we consider the case A << Tk- As we lower the energy scale to A, the operator /o or / is frozen, but /i is 
proportional to A _JV °/ 4 : 

/ 1|t = doA-^V + .... 

The operator g represents the single electron level at zero energy (the number of electron levels is odd, since fo is 
frozen) . 

Now the effective Hamiltonian (at the energy scale A) can be written as 

H eff = -JS-a(0) + Kf+fr + w(f + g + h.c.) . (6) 

The operators / and g arise when we express fo and /i in terms of the lowest single electron/hole levels of the 
conduction electron Hamiltonian. In the effective Hamiltonian given above J and K are renormalized coupling 
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constants; they increase in magnitude as the high-energy degrees of freedom are integrated out but they saturate at 
the value attained at an energy scale of Tk and are not further altered (since / is frozen); however, the coupling 
constant w will continue to scale as w oc hr N °/ i when Nq increases (or as A decreases), we expect w oc y/A. Since w 
should be of the order of Tk when A = Tk, we can re- write w = olTk\J A/Tr-). Note that in writing down the above 
Hamiltonian we have neglected all irrelevant terms, the inclusion of which will not change the results qualitatively. 

We want to investigate the nature of the ground states of the above Hamiltonian for the cases K = and K ^ 0, 
by diagonalizing the Hamiltonian. This is mildly tedious but can be carried out in a straightforward fashion. The 
main results are as follows: The ground state is always a singlet (when K > 0, the ground state is in the subspace 
(Q = —1, 5 = 0) (for K < 0, it is in the subspace (1, 0)). The first excited state is in the subspace (0,1/2) and has a 
gap relative to the ground state proportional to A. For K = 0, the ground state is in the subspace (0, 1/2), which is 
a doublet with the energy gap to the first excited state proportional to A 2 . These results are in agreement with our 
numerical RG computations. For the benefit of the reader a derivation of these results is presented below. 



A. Diagonalizing the effective Hamiltonian 



We diagonalize H e ff in Equation (|fj) in two steps. Diagonalizing the first two terms of the Hamiltonian H e ff in 
the subspace of the / states gives rise to four eigenstates given below: 

• State A (-1,1/2): E = 

• State B (0,0): E = §J + K 

• State C (0,1): E = -\J + K 

• State D (1,1/2): E = 2K 

Here the numbers in the parentheses denote the charge and spin of the energy states. 

Now we add the g states. The Hamiltonian can be written in the basis consisting of A, B, C, D and g states using 
a procedure similar to what was employed in the iteration-scheme of Wilson's RG iteration (see for example, Eqn. 
(B2) in Appendix B of the paper by Krishna-murthy et alEl). Let Al, A2, A3, A4 denote the basis states obtained 
by combining A with zero, one, and two g states, etc. The Hamiltonian matrix in each charge-spin subspace can be 
written down as 

• State Al (2,1/2): H 2 ,i/ 2 = 

• State A3+B1 (-1,0): 



State A2+C1 (-1,1): 



H-i,i = 



w -\J + K 



State A4+B2+C3+D1 (0, 1/2): 



H, 



0,1/2 



f 73 



V 





w 



-\J + K 



V2 



2 

2K 



State C2 (0,3/2) 



H 



3/2 



K 



G 



• State B4+D3 (1,0): 

Hi.o 

• State C4+D2 (1,1): 

#1,0 

• State D4 (2,1/2): 



\J + K -w 
-w 2K 



■\J + K -w 
-w 2K 



#2,1/2 = 2K 

Whether the ground state is a singlet or doublet depends on the relative energies of the lowest energy levels in 
subspaces (—1,0), (1,0), and (0, 1/2). If the lowest energy level in the subspaces (—1,0) and (1,0) is lower than the 
lowest energy level in the subspace (0,1/2), then we have a singlet (T\ will approach zero); otherwise, we have a 
doublet and T\ approaches 1/4. 

Let us first consider the case K ^ 0. We perform a second-order perturbation calculation of the energy of the 
eigenstate with the eigenvalue near | J + K: 

• For the subspace (—1,0), we have 



p + K 



For the subspace (1,0), we have 



For the subspace (0, 1/2), we have 



E « ^J + K 



\J-K 



3 T T ^ 1 w 2 1 w 2 
E « - J + K + --5- + - 



2\J + K 2\J-K 

It is clear that the ground state is always a singlet: when K > 0, the ground state is in the subspace (—1,0), 
whereas for K < 0, the ground state is in the subspace (1,0)). The energy level of the first excited state (in subspace 
(0, 1/2)) relative to the ground state is (assuming Kq > 0) 



w 2 1 w 2 



2\J-K 2\J + K" 



which is proportional to A. The energy level of the second excited state (in subspace (Q = 1, S = 0)) is E2 ~ 2E\ 
(this result was also found in our numerical results for the energy levels). 

How about Ko = 0? The issue cannot be resolved at the level of second-order perturbation theory. A fourth-order 
perturbation calculation for the lowest energy in the subspace (0, 1/2) yields 

^°^4 J+ 3T + 27 J( J ) ' 
For the subspaces (—1,0) and (1,0), the lowest energy is given by 



C 3 7 // 3 t\ 2 , 2 3 4w 2 1 W 4 

Eo = A J -fA J) +W 2 *4 J+ 3T-2 J( 7 ) - 

It is clear that the ground state is in the subspace (0,1/2), which is a doublet. This agrees with our numerical 
results. The energy gap of the first excited state is proportional to w 4 or A 2 . 



7 



B. Response Function 



Let us consider the calculation of T « S Z ;S Z » in the ground state when A << Tk- By definition 

<< S z ; S z »= / < S z (t)S z > dr. 
Jo 

Close to zero temperature, we can write 

^ <? • <? ^- \<G\S z \I>\ 2 (l-cM-m-E G ))) 

« s z ,s z »- ^ — — , 

where \I > represents many-body states of the system and \G > denotes the ground state. For temperatures much 
smaller than the energy gap between the first excited state and the ground state, we have (separating out the 
contribution of the ground state from the summation) 

«S Z ;S Z »=(3\<G\S Z \G>\ 2 + J2 1< e^e >12 - 

\I>*\G> ' G 

Since the second term in the above expression is finite, we obtain in the limit as T — > 0, 

T«S Z ;S Z »=\<G\S Z \G>\ 2 . 

For the case Kq ^ 0, it is easy to verify that T << S z ; S z >> is zero. 
For the case that Kq = 0, we find that 

4w 2 

<G \S Z \G>=±—. 

Thus T << S z ; S z >> in this case is proportional to u> 4 or A 2 in agreement with the numerical results. 



C. Anderson Model 



We now discuss the Anderson Model briefly since the results are similar to those of the Kondo problem discussed 
above. We consider the limit that U is very large and T is very small. The effective Hamiltonian is of the form 

H ef f = e d d+d^ + V{d+f Qlx + h.c.) + A( a + &Al - b+b^). 

Here €d is the effective impurity level at the energy scale A and V is the effective coupling to the conduction electron 
states. The impurity level cannot be doubly occupied: < 1. For the case that — > A, then the local moment 
regime will be reached, and the effective Hamiltonian can be converted to the Kondo Hamiltonian; this has been 
discussed above. Here we focus on the case that — a « A. 

Consider the case when V is very small; to leading order, the ground state depends on the sign of e<j- If £d > 0, then 
the ground state corresponds to two electrons occupying the conduction electron level at —A, and it is a singlet (this 
situation arises, for example, when the initial is greater than zero; this has been checked by the nonperturbative 
RG calculation). If < 0, then the ground state corresponds to two electrons occupying the conduction electron 
level at —A and one electron occupying the impurity level. So the ground state is a doublet. 

From the above analysis, it is clear that the reason for the ground state being a doublet for the case — > A (when 
the local moment regime is reached) and — << A (when only the mixed- valent regime is reached) are different. 
In the first case, the local moment regime is reached, and as the temperature is lowered, the moment begins to be 
quenched due to large effective \J\, but as the temperature is further lowered, one can see the small splitting of the 
singlet state to the doublet state due to the finite gap energy. In the second case, the local moment is not formed at 
the energy scale A. But as the temperature is lowered to the energy scale of — e^, charge fluctuations are suppressed 
and they eventually cease to exist, and the system becomes a doublet. 
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V. CONCLUSIONS 



We have performed a Wilson Renormalization group calculation of the Kondo and Anderson models with a gap 
in the conduction electron density of states. The impurity susceptibility correlation functions, and a zero-frequency 
response function have been calculated as functions of temperatures in various regimes. Our calculations confirm 
earlier results on the qualitative differences in the low-temperature behaviors between the cases with and without 
particle-hole symmetry when the gap is much smaller than the Kondo temperature. We have shown that the numerical 
results at low temperatures can be understood in terms of simple low temperature effective Hamiltonians. 
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FIG. 1. Tximp plotted as a function of T for the Kondo problem. The gap energy is A = 1.22 x 10~ 4 . (a) The potential 
scattering is absent (K = 0.0). The values of the coupling J used are —0.1, —0.2, —0.3, —0.4, and —1.0. Note that as T — > 0, 
T\ approaches 1/4. (b) Particle-hole symmetry breaking is present (K — 0.1). Note that for J = —0.4, T\ goes to zero, while 
for J = —0.2, it approaches 1/4. 

FIG. 2. Tximp plotted as a function of T for the Anderson model, (a) Symmetric Anderson model with U = 0.1, e d — —0.05, 
and T = 0.006. Note that as T — ► 0, T\ approaches 1/4. (b) Asymmetric Anderson model with U = 0.1, en = —0.001, and 
r = 0.00015. Note that as T — > 0, T\ approaches 1/4 if A >> T K and if A << T K - The values for the gap energy A used 
in the calculations are shown in the legends of the figures. 

FIG. 3. The zero-frequency response function T « S Z ;S Z » plotted as a function of T for the asymmetric Anderson 
model with U = 0.1, e d = —0.001, and T = 0.00015. Note the qualitative differences between T << S z ; S z » and T\ at low 
temperatures. 

FIG. 4. The local correlation functions: (a) < nj > and (b) < S ■ <r(0) >, plotted as a function of T for the asymmetric 
Anderson model with U = 0.1, e d — —0.001, and T = 0.00015. Note that when A << Tk, the correlation functions approach 
constant values independent of A as T — > 0. 

FIG. 5. The impurity susceptibility (Fig. 5(a)), < n d > (Fig. 5(b)), and T << S Z ;S Z >> (Fig. 5(c)), plotted as a function 
of T for the asymmetric Anderson model with U = 0.1, e d = —0.025, and T = 0.01. 
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TABLE I. The values of T « S z ; S z » at zero temperature for a range of values of the band gap. The numbers enclosed 
in the parentheses are the total charge Q and spin S of the ground state. 



A 


Jo = -0.2, K = 0.0 


Jo = —0.2, K — 0.0 


1.88 x l(T b 


0.0398 (0, 1/2) 


0.0490 (0, 1/2) 


6.27 x lCT 7 


0.00797 (0, 1/2) 


0.0110 (0, 1/2) 


2.09 x 10"' 


0.00101 (0, 1/2) 


< 1.2 x 10~ b (-1, 0) 


6.96 x 10~ s 


0.000115 (0, 1/2) 


< 4.0 x 10 _b (-1, 0) 
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